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Abstract 

We explore the viability of baryogenesis from light scalar decays after the 
electroweak phase transition. A minimal model of this kind is constructed 
with new CP violating interactions involving a heavy fourth family. The 
departure from thermal equilbrium must come from topological defects like 
cosmic strings, and we show that almost any mechanism for producing the 
cosmic strings at the electroweak scale results in a viable theory. Baryogenesis 
occurs in the fourth generation but the baryon number is later transported to 
the visible generations. This mechanism of indirect baryogenesis allows us to 
satisfy experimental limits on the proton lifetime while still having perturba- 
tive baryon number violation at low energies. The fourth family has very small 
mixing angles which opens the possibility of distinct observable signatures in 
collider experiments. 
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1 Introduction 



The experimental bounds on the ratio of baryon excess to the entropy of the universe are 

r/ = (2 - 8) X 10-^\ (1) 

Baryogenesis is an attractive explanation of the observed fact that baryons are more abun- 
dant than anti-baryons in the universe. The conditions necessary for baryogenesis were 
spelled out by Sakharov nearly three decades ago 0. The three conditions are: (i) exis- 
tence of baryon number violating interactions; (ii) C and CP violating processes; and (iii) 
a departure from thermal equilibrium. The earliest models of baryogenesis were based on 
baryon number and CP violating processes of GUT theories. The necessary departure 
from thermal equilibrium was achieved by having superheavy bosons decay by slow inter- 
actions that make them overabundant (in comparison to their thermal distribution) in a 
rapidly expanding universe |P, ^. 

It was however realized subsequently that anomalous baryon number violation in 
the electroweak theory itself could wipe out any baryon asymmetry formed at the 
GUT scale unless the density of baryons minus leptons {B — L) is non zero. Another 
difficulty with GUT scale baryogenesis is inflation. Inflation is needed to get rid of heavy 
monopoles formed during GUT scale symmetry breakings, but it also inflates away any 
baryons produced at that scale. 

Since then several other mechanisms for baryogenesis have been proposed that produce 
baryons at or after the electroweak phase transition. The most notable is the mechanism 
of electroweak baryogenesis This mechanism, however, requires a sufficiently strong 
first order phase transition 0, |^ . At present the question of the order of the electroweak 
phase transition remains unanswered. If the electroweak phase transition were weakly 
first order or second order, then alternative mechanisms of baryogenesis at or below the 
electroweak scale would become very attractive. 

Such models can be classified broadly by asking the two fundamental questions: 

(i) What is the source of B violation? 

(ii) What is the reason for departure from thermal equilibrium? 

The question oi CP violation is not included among the above two criteria. Models 
of baryogenesis must go beyond the standard model to incorporate new CP violating 
interactions. However, we do not see any obvious way of classifying the new CP violating 



1 



sectors. 

The usual source of B violation is one of the following two: 

A) The non-perturbative B violation in the standard model. 

B) Perturbative B violation in an extended standard model. 

We would broadly classify the reason for departure from thermal equilibrium into:0 

a) Out of equilibrium decay of massive excitations. 

b) The expanding or collapsing wall. (The term "wall" refers to a sharp change in the 
value of the Higgs field's vev.) 

The mechanism of electroweak baryogenesis falls into the classes A) and b) respec- 
tively. The "wall" is the expanding wall of a bubble of true vacuum at the electroweak 
phase transition. If we are looking for alternatives to electroweak baryogenesis but want 
baryogenesis to occur after the electroweak phase transition, then we either need to give 
up one or both of A) and b) or look for a model that falls in the classes A) and b) but 
does not require a first order phase transition. 

There are in fact two mechanisms in the literature that are similar to electroweak 
baryogenesis in that they fall in the classes A) and b) but the "wall" is provided not by 
an expanding bubble of true vacuum but by a collapsing cosmic string. In one of them , 
anomalous baryon number violating processes take place inside collapsing cosmic strings. 
In another case [rU|, electroweak strings (or Z strings [^) with magnetic Z flux are 
needed. Both these models need new CP violating sectors. The first one also requires 
that sphaleron effects be appreciable inside the "core" of the cosmic strings. All models 
of this kind are therefore sensitive to the structure and properties of the strings. In some 
cases they may require a light Higgs which implies a first order phase transition so the 
question of finding an alternative to electroweak baryogenesis remains open. The second 
model needs metastable Z strings, for which no viable extension of the standard model 
exists so far. 

The limitations of the above models lead us to explore other alternatives to electroweak 
baryogenesis where we give up one or both of A) and b). The GUT scale baryogenesis from 
heavy boson decays in fact falls under the classes B) and a). These models automatically 
have the virtue of being insensitive to the order of the electroweak phase transition. 
The question arises, whether there are viable models falling under these classes that can 
produce baryons below the scale of the electroweak phase transition. 
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Other mechanisms have been suggested. For example see ref. [12| 
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Seemingly, there are two hurdles to having this mechanism work at such low energy 
scales. The first is the constraint from proton decay. The stability of the proton implies 
that baryon number violating interactions must couple to the first generation quarks with 
a very small coupling. It seems that to make this ratio small the baryon number violating 
interactions must involve superheavy particles (mass > lO^^GeV). The second obstacle is 
that for a massive excitation to decay out of equilibrium, its decay (and annihilation) rates 
should be smaller than the expansion rate of the universe. With SU{3)c x SU{2)w x 
U{1)y couplings or with Yukawa couplings of the order of 10^^ — 10~^, the decay (or 
annihilation) rates are usually greater than the Hubble expansion rate unless the universe 
is at temperatures as high as 10^'' GeV. In ref. [0 these two hurdles were overcome 
by a model where proton decay was forbidden by lepton number conservation and some 
heavy excitations were required to be SU{3)c x SU{2)iy x f/(l)y singlets. The model 
needs colored scalar s as well as a pair of massive Major ana fermions which have no gauge 
interactions. In this paper we seek an alternative particle decay mechanism with, what 
we believe, a simpler spectrum that may still have distinct experimental signatures. 

Our motivation is in fact twofold. Firstly we would like to incorporate the merits of 
boson decay models in a model that can be completely described by an effective theory 
at the electroweak scale. By bringing down the scale of B and CP violation to the elec- 
troweak scale one improves the testability of the theory compared to GUT scale models. 
We would also like to restrict the fermion content to simple sequential families and would 
not require any of them to be gauge singlets. Naturally the model should be viable even 
if the electroweak phase transition is second order. 

Secondly, extensions of the standard model, such as models of top-color assisted tech- 
nicolor [0 , suggest the possibility of having symmetries under which quarks of different 
families transform differently. These symmetries must be broken above the electroweak 
scale to permit quark mixing at low energies. However the existence of these symmetries 
opens up the possibility of having very small quark mixing angles. Thus the problem of 
proton decay that must be addressed in models of baryogenesis through perturbative B 
violation may find a new solution through small mixing angles (10~^ — 10~^) between new 
heavy quarks that have B violating interactions and the light quarks that constitute the 
proton. If a viable model of baryogenesis makes use of this mecahnism, the extra fermions 
needed in the CP violating sector would simply consist of copies of the observed fermions, 
yet can have distinct signatures in collider experiments which would point strongly to an 
underlying theory with perturbative baryon number violation at low energies. 
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In this paper we show that the idea of small mixing angles described above does indeed 
provide a natural answer to the needs of any model of baryogenesis based on light boson 
decays. In our model, although baryogenesis happens through the decay of scalar bosons, 
the departure from thermal equilibrium can be obtained only by having topological defects 
like cosmic strings. This is a consequence of making the theory insensitive to the order 
of the electroweak phase transition . However, we show that once the simplest B and 
CP violating sector incorporating the above ideas is constructed, the cosmic strings can 
be obtained as an inevitable bonus. In addition, the mechanism is not sensitive to the 
structure and properties of the cosmic strings. 

The paper is organized as follows. In section 2 we lay out the basic picture of how 
the model works. In particular we point out the various parts of the mechanism that 
one must check carefully to compute the baryon asymmetry generated in the end. The 
phenomenological viability of the model is then shown part by part in the following 
sections. 

In section 3 we describe the CP and B violating interactions and review the boson 
decay mechanism of 0] that forms the core of the present mechanism. We then consider 
phenomenological constraints on the mixing angles of the model from proton decay ex- 
periments. Finally we comment on the potentially observable experimental signatures 
predicted by the model. 

In section 4 we present an estimate of the number density of the scalar bosons generated 
by the decay of cosmic strings in this model. 

In section 5 we consider the Boltzman's equations for the evolution of baryon number 
after the electro- weak phase transition and show that for a particular range of parameters 
the baryon asymmetry generated immediately after the phase transition can survive till 
the present time. 

2 Basic Mechanism 

Our starting point is the boson decay mechanism. Simplest models of baryogenesis need 
at least two bosons with B and CP violating interactions. In the next section we will 
give a brief review of this mechanism. As mentioned earlier, for successful baryogenesis 
there must be departure from thermal equilibrium. However our primary motivation is 
to avoid having superheavy scalar bosons. In the usual picture of cosmological evolution, 
TeV scale excitations do not go out of thermal equlibrium unless they are practically 
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non-interacting. Therefore we are led to consider the following scenario. Suppose cosmic 
strings are formed at or slightly above the electroweak phase transition. Some of these 
strings will be in the form of loops which will subsequently decay into particles. If a large 
number of these loops decay into scalar bosons that are sufficiently heavy, then there may 
be an overabundance of the scalar bosons. The decay of the overabundant scalar bosons 
can then be the reason for the necessary departure from thermal equilibrium. 

We would like to construct the minimal B violating sector of scalar bosons of the 
above kind. As shown in ref. ^ there will be at least two of these bosons. If they have 
direct B violating couplings to light quarks then we have a large width for proton decay 
unless the couplings are very small. But if the couplings are very small then we do not 
seem get sufficient B violation and there is no baryogenesis. 

A possible solution of the above problem is to have a fourth generation of quarks. 
The fourth generation can also be used to include new CP violation which is a necessary 
ingredient of baryogenesis. The scalar bosons can decay into the new quarks through B 
violating processes leading to baryogenesis. However the fourth generation quarks must 
mix with the lighter generations for baryogenesis to occur in the visible generations. In 
the next section we show that it is possible to have small mixing angles (< 10~^) between 
the B violating fourth generation and the other three generations such that the width for 
proton decay is within experimental bounds. 

Thus the basic picture is the following. Immediately after the electroweak phase 
transition there must be a network of decaying string loops that produce a thermally 
overabundant quantity of scalar bosons. These scalar bosons decay into a fourth family 
of quarks and leptons while generating a baryon asymmetry. Finally the fermions of the 
fourth family must decay into fermions of the other families. This is the minimal structure 
for baryogenesis from light scalar decays after the electroweak phase transition. 

There are several cosmological and phenomenological constraints that the model must 
satisfy. Here we list the main constraints that will be shown to be satisfied by the model 
in the next sections. 

(i) All couplings must be natural. We do not address the hierarchy problem due to 
scalars that also exists in the standard model. 

(ii) Fermions of the fourth family must have experimentally viable masses. 

(iii) There should not be a large width for proton decay. This implies small mixing 
angles between the fourth and the other families. 
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(iv) Cosmic strings should be naturally incorporated in the theory. 

(v) All heavy particles (the scalar bosons and the members of the fourth family) must 
have large decay widths so that they decay within a cosmologically acceptable period. A 
long lived heavy particle may cause the problem of the overdosing of the universe. In 
particular, the fourth generation leptons must mix with the other three generations in 
order to decay. The lepton mixing must be large enough from this point of view. 

(vi) The baryon asymmetry, once created, should not suffer a wash out from inverse 
decays. 

Each of the above constraints is fairly restrictive. Indeed it is not obvious that (iii) and 
(v) can be satisfied simultaneously. However, as we show now, the minimal model built 
with the motivations that we have mentioned earlier satisfies all the above conditions for 
viability. The rest of the paper is devoted to the viability proof of this mechanism. 

3 B Violation and a Fourth Generation 

3.1 B and CP Violation 

The minimal hght scalar boson decay model must have: 

(i) two scalar bosons Xi and X2; 

(ii) a fourth family of quarks and leptons (1/4, 64)^, 1/411, ^ar, {t', b')L, t'u, b'j^. 

Note that we have a right handed neutrino in the fourth family which is required by 
experimental bounds on the number of hght neutrino flavors (2.99 ±0.04) in the standard 
model. We can have all the necessary B and CP violation with the above particle content 
if Xi and X2 are SU{2)w singlets and SU{3)c triplets and have hypercharge Y — —2/3. 
In addition to the coupling of all the fourth generation quarks and leptons to the standard 
Higgs, we can then have the following new Yukawa couplings: 

Ly = ...giit'jib'jiXi + gnh'RV^Xi + gi^t'j^e^Xi + giJALqlEXi + gm^LqiLXi 

+h.c + .... (2) 

with / = 1, 2 and representing the new quark doublet. 

With these interactions the action has a non anomalous Yl\=i{.B — L)j x [B — L)^ 
symmetry, where Yll=\{B — L)i is the difference of the baryon and lepton numbers in the 
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first three generations and [B — L)^ is the difference of the baryon and lepton numbers 
in the fourth generation. 

The CP violation in Ly can be communicated to the ordinary quarks and leptons if 
the quarks and leptons of the fourth generation mix with the other generations. Note that 
lepton mixing is allowed since the neutrino in the fourth generation is massive. When 
these mixings are present, the quarks and leptons of the lighter generations must couple 
directly to the colored scalars Xi and X2 and the global non-anomalous symmetry is 
reduced to {B — L). We require that terms in the action that violate {B — L)^ by n units 
are suppressed by small couplings and mixings of the order of s" where s is a small number 
~ (10^"^ — 10"^). One way to motivate this suppression is to think of this action to be a 
low energy limit of a theory that has a. {B — L)^ violating sector involving only massive 
fields. When these massive fields are integrated out, the resulting {B — L)^ violation is 
suppressed by powers of their mass. In Appendix A we describe an explicit way of realizing 
this suppression by considering a simple theory with gauged {B — L)i^ — J2i=i{B — L)i. 
When this symmetry is broken, a low energy theory similar to the one described above 
is obtained. Other allowed couplings are gauge invariant quartic couplings between the 
new scalars and the standard Higgs. The suppression rule described above should apply 
to all these couplings. 

The interactions (§) violate baryon and lepton numbers while preserving B — L. The 
interactions are in fact identical to the scalar boson interactions of ref. [^, where it was 
shown that the out of equilibrium decays of Xi and X2 produce a baryon asymmetry. 
The amount oi CP violation coming from the above terms was calculated in and we 
briefly recapitulate the main results. At tree level, decays of Xi and X2 produce no CP 
violation since the cross section for Xi — > baryons is exactly equal to the cross section 
for Xi —>■ antibaryons. The same statement applies to decays of X2 and X2. However at 
the one loop level there are several other processes contributing to the same decay modes. 
Consider for instance the decay Xi —>■ t'^ + b'j^ and Xi — > + z/4/j. The one loop diagrams 
contributing to the process have an internal X2 propagator (Figure 1). The interference 
of the tree order and one-loop diagrams produces a net CP violation for the decay of 
Xi through these channels. The baryon number produced per decay through the channel 



where / is the relevant loop integral. If all the fermions have similar masses, the loop 
integrals in all the channels are the same. Taking all channels into account one then has 



X, 



AB 



4Im(^ii^*i^i2^22)Im(/) 

9u9*n 



(3) 
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the expression 



AS = \g\Hm{I)e, 



(4) 



where g is a, typical Yukawa couphng and e is a phase angle characterising the average 
strength of the CP violation. The imaginary part of / is easily evaluated when all the 
fermions are massless: 



where p — mxi/mx2 is the ratio of the masses of the two colored scalars. Clearly the 
CP violation is zero if the two scalars have the same mass. For comparable but unequal 
masses Im(7) ~ (10~^ — 10~^). The value of Im(7) decreases if the fermions are not 
massless, but the order of magnitude estimate is unchanged. 

The value of e can be as large as 1 radian. Once a range of values for the masses 
and mixing angles are found, the allowed range of values for e is fixed by the baryon to 
entropy ratio generated in the theory. 

3.2 Masses, Mixing Angles and Proton Decay 

We shall see later that for baryogenesis to be successful it must be possible for the fourth 
generation quarks and leptons to decay into lighter quarks and leptons, while the reverse 
process must be prohibited. Therefore all fourth generation fermions need to be heavier 
than the Z, so that they can decay to a or a Z and a hght fermion. The mass difference 
between two members of an electroweak doublet can not be large by considerations of the 
p parameter. In this model, since all masses come from standard couplings to the Higgs, 
there is no obstruction to satisfying this criterion. Masses of the scalars Xi and X2 need to 
be slightly higher than the electroweak scale. A mass of a few TeV seems to be necessary 
for sufficient baryogenesis. 

The mixing of the fourth generation quarks and leptons to members of the other three 
generations provides a way for the heavy quarks and leptons to decay. As we shall see 
later, with too small mixings baryogenesis never occurs in the visible sector. On the 
other hand if the mixings are too large, the decay width of the proton increases beyond 
experimental bounds. We show that there is an allowed region in the space of the mixing 
angles where baryogenesis is achieved while having acceptably small decay width for the 
proton. 





8 



In GUT theories where quarks of the first generation have baryon number violating 
couphngs, the proton decays by processes shown in Figure 2. To meet the experimental 
limit on the lifetime of proton (> lO^^years), the mass of the internally propagating GUT 
boson should be > 10^^ GeV. 

In the present model B violating couplings involving the first three generations are 
down by a suppression factor. For instance, if the mixing angle between a fourth gen- 
eration quark and a first generation quark is ©94^1, a typical baryon number violating 
coupling involving first generation quarks will be down by two powers of this mixing: 

Ly = .... + gnOlg^URdRXi + .... (6) 

We retain the coupling gn to exhibit that the couplings now are smaller in comparison 
to the similar couplings involving only fourth generation quarks. The processes leading 
to proton decay still look like the one sketched in Figure 2, except the baryon number 
violating vertices now have very small coupling constants. The amplitude for the process 
is 

^~/©l..©W3-V, (7) 

where Qq^q^ is a typical mixing angle in the quark sector (between the first and fourth 
generations), Qm-j, is the largest mixing angle in the lepton sector and is the product 
of the typical coupling constants at the two baryon number violating vertices without the 
suppression factors. 

Using mx ~ 10^ GeV, we find the following inequality for the mixing angles if exper- 
imental bounds on proton decay are to be satisfied: 

l/e^..,ew,|<lo-^^ (8) 

In principle it is possible to have greater quark mixings if the lepton mixing is smaller. 
Also, ©5453 and ©5452 are likely to be greater than Qq^q^ by one and two orders of magnitude 
respectively. 

The spectrum of masses and the mixing angles allowed for the fourth family makes 
for interesting and distinct experimental signatures. The presence of a sequential fourth 
generation is not ruled out by experiments. The D0 collaboration puts a limit of rUg^ > 
131 GeV from the charged current (CC) decay modes of the fourth generation quarks t' 
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and h' (t'^h^W^h'^t^ W\ However if the h' is lighter than the top quark but 
heavier than the Z, then its dominant decay mode is the flavor changing neutral current 
(FCNC) mode 6' ^ 6 + Z A search for h' through this mode is currently feasible and 
one expects some experimental results in the near future (see ref . |T7[] ) . 

The z/4 is also potentially observable through the tri-lepton decay mode — > Zz//i. 
The LEP I bound on the heavy neutrino mass remains niy^ > 46 GeV. Some restrictions 
on the mixing angles of the neutrino have been placed by the D0 collaboration [^] which 
are consistent with the mixing angles allowed in the present model. As we shall show 
later, the mixing angles between the fourth and the third generation quarks and leptons 
can be in the range 10^^^ < 65453 ~ '^i^i-^ < 10~^. Very small mixing angles will make 
the lightest fourth generation quark effectively stable inside the detector and may lead to 
peculiar signatures. Even if the mixing angles were measured to be closer to the upper 
bound of the above range, the obvious natural explanation for their smallness would be 
the existence of extra symmetries at scales higher than 1 TeV which forbid the mixing 
between the fourth generation and the other generations. Then, with the standard gauge 
and Higgs couplings, there is an extra symmetry in the form of B4 — L4 above the TeV 
scale. For the small quark mixings in the fourth generation to exist, this symmetry must 
be broken at the TeV scale and a likely result would be perturbative violation of baryon 
number in the fourth generation. Therefore experimental signatures that are consistent 
with the mixing angles predicted in this model point very strongly to a mechanism of 
baryogenesis through perturbative B violation. 

4 Cosmic Strings and Baryogenesis 

With the B violation and CP violation in place, all we need to produce baryons is an 
out-of-equilibrium decay of the scalar bosons. To achieve this there has to be a mechanism 
for making the scalar bosons overabundant immediately after the electroweak scale. This 
can happen through the formation and decay of cosmic strings. The cosmic strings must 
form close to the electroweak phase transition. Strings formed much earlier will have a 
distribution with a correlation length that is too large and their number density will be 
too small to generate the baryon excess we see today. Since the symmetry breaking in 
the standard model does not produce any cosmic strings, extra broken gauge symmetries 
must exist. However, as described in Appendix A, the broken symmetry can be the 
approximately conserved {B — L)4 — J2i=i{B — L)j. If this is a gauged symmetry that 
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is spontaneously broken at the electroweak scale, the smallness of the mixing angles and 
the presence of the cosmic strings can be explained simultaneously. 

Moreover, besides being candidates for seeding galaxies, cosmic strings occur in many, 
independently motivated extensions of the standard model. For instance the Aspon model 
has an extra U{1) that breaks close to the electroweak scale. The motivation there is to 
provide an explanation for the small value of the vacuum angle 6 in QCD. Supersymmetric 
extensions of the standard model have been proposed which attempt to resolve the fj, 
problem of the minimal supersymmetric standard model (MSSM) and the cosmological 



solar neutrino problem and which have strings |^0|- Both these models were considered 



recently |21] in the context of electroweak baryogenesis from cosmic strings. Top color 
models |14| are another class of models where a U{1) gauge group is broken close to 
the electroweak scale. In our case, almost regardless of motivation, any extension of the 
standard model that produces strings at the electroweak scale will work. The cosmic 
strings need not have a special structure or satisfy any particular requirements. The 
effective scalar whose vacuum expectation value (vev) causes the formation of the cosmic 
strings will naturally have quartic couplings with the bosons Xi and X2. This ensures 
that Xi and X2 will be produced from the decay of cosmic string loops. 

The scenario for baryogenesis is now very similar to where the authors considered 
emissions of heavy particles from collapsing cosmic strings. If the heavy particles are 
produced at a scale which is sufficiently small compared to their mass, they may become 
overabundant and through CP violating decays generate baryon number. In the present 
model we focus on the overabundant production of Xi and X2 particles from strings. 

Immediately after the electroweak phase transition, the space will be filled with a 
criss-crossing of string network that looks like a random walk in three dimensions. The 
initial correlation length is ipt^ ~ Numerical simulations indicate |23] that a large 
fraction (~ 80%) of the total string length resides in the infinite strings. The rest is in 
the form of loops which have a scale invariant distribution 

— = R-^ (9) 
dR ' 

where R is the characteristic size of the loops. The initial network has loops which decay 
rapidly. The infinite segments also generate more loops by frequent intercommutations. 
The net result is that the correlation length increases with time and the string network 
enters a scaling solution when the correlation length equals the horizon size p^ . 

In the literature the period from the time of the phase transition (tc) to the time 
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t* = {Gfi)~^tc is called the friction dominated period [G is the Newton's constant 
and fi is the mass per unit length of the strings). Loops produced at time t with tc <t <t* 
immediately shrink to the radius 

r;(t) = G/imJ/V/^ (10) 

where mpi is the Planck mass (~ lO^^GeV). Below this scale friction effects are sub- 
dominant and the loops shrink predominantly by processes like gravitational radiation 
and cusp annihilation. The shrinkage rate from gravitational radiation is [EH 



f = -7G.. (11) 

where 7 is numerically determined to be about 10. 

The time for a loop to shrink to a size which is of the order of its thickness due to 
purely gravitational effects is to = {'yGfi)~^R. In the present case tc ~ lO^^R ~ lO^^s 
even for the smallest loops {R ^ ip). This is a very long time compared to the cosmological 
times of interest and gravitational effects are therefore completely negligible in our model. 

Cusp annihilations, on the other hand, can occur at a much faster rate. The rate of 
shrinkage in this case can be modelled by [^] 

dR 7p 



dt {Rr]y/^ ■ 



(12) 



The corresponding decay time is tcusp = ^^^"^^ . As shown in this is less than one 



expansion time if the time of formation t < {-Q^Y^'^tc- For Gfi ~ 10~^^ it is safe to assume 
that this condition is satisfied for a long time. Thus in our model, the loops formed at 
time t immediately shrink to R = rf{t) after which they shrink by cusp annihilation to 
R ~ within one expansion time. The lifetime of a loop formed at time t is then 

r(t) ~ t. 

When the string loops have shrunk to a radius that is comparable to their thickness 
^-1/2 (where fi is the mass per unit length of the string), nonlinearities in the scalar field 
potential will cause the entire loop to decay into elementary particles. It is from this final 
burst process that we can expect the heavy scalars Xi and X2 to be produced. Since 
/i ~ m?x^^m?x^\ the number of Xi and X2 particles that we get from each loop is ~ 1. 

The baryon number produced by string decays can now be evaluated by computing 
the number of string decays from the time t^. of the phase transition up to the present 
time. There are two kinds of loops to be taken care of; those formed at the time of phase 
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transition {t = t^) and those formed after the phase transition (t > t^). Also, the baryon 
number produced from the decay of Xi and X2 bosons will be washed out unless baryon 
number violating processes are effectively frozen after the Xi and X2 decays. To compute 
the remnant baryon asymmetry we therefore need to evaluate the number of string decays 
that take place after t = tf, where t/ is the freeze out time for baryon number violating 
interactions. Two cases, thus, arise: 

1) tf^t. 

The number of decaying loops is the sum of loops produced at tc and after t^.. The 
number density of heavy bosons produced from loops formed at t^ is obtained by inte- 
grating (1): 



where the last factor takes care of the dilution of the number density due to the expansion 
of the universe. More cosmic string loops are formed after the phase transition as the 
coherence length increases with time and loops are chopped off from 'infinite' sections of 
the string network. The loop production rate is related to the rate at which the coherence 
length increases by |22| 

^ = (U) 
dt ^Ij^ dt ' ^ ^ 

where u is a. constant of the order of unity. Integrating (14) we have 



".>.. = A' ^P^M'/''~iV7^«"... (15) 



, —r[ — ] # ~ N ^ , 

In order to have a large n^^ we need to minimize V'tc- Using ■i/'i^ ~ ^ for strings produced 
at the electroweak scale, we get the ratio of the total number density of X^ and X^ to the 
entropy density to be 

C^max ~ — ^ — , (16) 

s 

where s is the entropy density of the universe and g^, ~ 100 is the number of effectively 
massless degrees of freedom at the electroweak scale. 

2) tf > t. 

Now the loops produced at t = decay within one expansion time and the resulting 
baryon number is washed out. Thus we need consider only the contribution from the 
loops decaying after tf: 
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In the friction dominated period i/j t2 therefore we have, u = cUmaxlfy)^- Thus 
unless tf = tc, there is some damping in the production of baryons. In the next section 
we show that it is possible to have tf ^ in our model. 



5 Approach to Equilibrium 

The baryon asymmetry generated from Xi and X2 decays will evolve according to the 
Boltzman's equations. A large number of particles and interactions are relevant. For 
each species of particles one gets a Boltzman's equation. The equations are coupled 
integro-differential equations and for an accurate estimate of the baryon number, one 
must integrate them numerically. Useful analytical estimates can, however, be made by 
making simplifying approximations that reduce the number of degrees of freedom (and 
hence the number of Boltzman's equations) to a few. 

The greatest simplification results from considering a single scalar boson X instead of 
Xi and X2. In the following we consider the most dominant interactions of this scalar 
boson. The relevant processes are: 

1. XX qq,GGJl etc. 

2. X ^4^4 

3. X -^qth 

4. 54^4 giZ4 

5. q^W 

6. h I3W 

X is the generic colored scalar; q and / refer to generic quarks of any generation while 
g4 and I4 refer to a quark and a lepton of the fourth generation. The W in interactions 5 
and 6 represents the W bosons of the weak interactions. We denote gluons by G. 

The processes 1 are dominant annihilation channels for the colored scalars. The pro- 
cesses 2,3,4 are baryon number violating processes. The processes 2,3 generate the baryon 
number while their inverse processes and processes 4 can wash out the produced baryon 
number. The processes 5 and 6 transport baryon and lepton number from the fourth 
generation to the third generation. We have omitted processes of the kind ^4 — > qlliqi 
to simplify the Boltzman's equations. These processes are prohibited if the quarks and 
leptons of the fourth generation have nearly equal masses. Inclusion of these processes 
does not change our main results in a significant way. The processes of baryogenesis 
and 'freeze out' can now be treated as two distinct stages in the evolution of the baryon 
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number. 

5.1 X production and decay 

We will assume that immediately after the phase transition the decays of cosmic string 
loops raise nx,the number density of X, to about T^. (Later we show that nx > 10~^Tj^ 
may be sufficient for baryogenesis) . Since Mx > Tc, the Xs are overabundant and their 
number will decrease rapidly through decays and annihilations. 

The Boltzman's equation for X is [|I| 



dn 



X 



dt 



+ 3Hnx 



DP 



f^\M{X^ij)\^~Uf,\M\{ij^X) 



fxfx\M{XX ^ - hf,\M{ij ^ XX)f 



;i8) 



where DP. 



Wa2...M,b2- = Hi / (2^)32^, / Wf2E^('^^) ^ ~ ^jPh,) is the phase space 

volume element, H is the Hubble constant, is the phase space density of species i and 
\M{ij.. ..kl)\'^ is the matrix element squared for the process ij.. —>■ ..kl. The matrix 
element is summed over initial and final state color, spin and fiavor degeneracies. The 
number density nx is the number density of all X particles regardless of color. 

When the final state particles are light, (|18]) reduces to |^ : 



dn 



X 



dt 



+ 3Hnx 



[nx 



n 



cqi 

xl 



(rx) 



nx 



{n^^f {a{XX^ij)) 



(19) 



where (Fx) is the total thermally averaged decay width of X averaged over initial color 
degeneracies, is the equilibrium density of X and {(j{XX — > ij)) is the thermally av- 
eraged cross section for XX annihilations (averaged over initial state color degenaracies) . 

The decay modes for X are: X — > q4,q4,qll4. The dominant annihilation channels are: 
XX qq, GG. Since X is very heavy we can use the zero temperature decay width for 
Tx to good approximation. The same is true for the annihilation cross section with an 
appropriate value for the cm energy. In Appendix B we have computed these rates. Our 
results are: 



{a{XX-^tj)) 



ZTT 

^ (47raqcd)^ 



(20) 
(21) 
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where g is a typical coupling of X to ^4 and I4. 

Because X is overabundant, nx ^ Ux- The reduction in nx is initially dominated 
by annihilation processes. The decays overtake annihilations when the X number density 
reaches the critical value 

_9\g\' Ml 

nxcTit — o 77~ ~ ■ l^^J 



2 (47raqed)2 ■ 

From this point onwards, the annihilations are quenched out and most of the Xs decay 
through the CP violating processes producing baryons. 



5.2 Freeze out 

Once a large number of X decays have taken place, there is an excess of baryons (and lep- 
tons) over anti-baryons (anti-leptons) in the fourth generation. Baryon number violating 
processes like inverse decays (^4^4,54/4 X) or 2 — > 2 processes {qiq^ Qih) will tend 
to wash-out this excess. Decays to Ws and third generation quarks and leptons will also 
reduce the baryon excess in the fourth generation (although preserving the total baryon 
excess) . 

To see if a freeze out can occur, we look at the Boltzman's equation for ^4: 
drir, 



3 / DP„ 



+ 2 



+ 



CJ4IJ4,I34U 
DPx,qiqA 
DP— 

DPq^^q^W 



fqJq,\M'{qq - fq-JjlWml ' 

fx\M{X ^ 9494)1' - fqJqAM{qm ^ X)p 
fx\M(X ^ 94/4)1' - fqJu\M{q,h ^ X)\' 
fq,\M{q, ^ q.W)^' - fqJw\M{qsW ^ q,)\' 



qq)\' 



(23) 



I M' (94^4 



94/4)!^ and |M'(g4/4 



g'4g'4)P have primes on them to indicate that the 



matrix elements do not include s channel contributions in which the intermediate X is 
on shell (a physical particle), since these contributions have already been included in the 
decay and the reverse decay terms. The full matrix element (squared), |M(g4g4 ^'4/4)1^, 
is related to {M'^q^q^ ^'4/4)1^ by 



|M(g4g4 54/4)1' = \M'{q4q4 q^U) 



+ 



TT 



MxTx 



-6 



p ri) +pi(2) - mj, \M{q,q, ^ X)\'\M{X qlh 



(24) 



Following |2^ we will simplify (]23|) by parametrizing the CP violation of the system 
in the following manner. We define the matrix elements Mq and the numbers r] and fj by 



|M(X ^ 5454) 



|Mop(l + r/)/2 =|M(gigI^X)p 
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\M{X^q^k)\^= |Mop(l-r/)/2 =\M{q,k-^X)\\ 
|M(X^gIgl)|2= |MoP(l + r7)/2 = \M{q,q, ^ X)\\ 
\M(X^q,k)\'= \Mo\\l-r])/2 = \M{qik ^ X)\\ (25) 

We have used CPT and unitarity to relate the squared matrix elements. Note that 
all matrix elements are summed over initial and final state spin and color degeneracies. 
The CP violation parameters rj and r] are related to AB by the relation 

AS = (r/-r/)/4. (26) 
Since all the quarks and leptons are in thermal equlibrium we have 
/,.(p) = e-^/^+-/^-e-^/^(l + ^) , 
Up) = e-^/^--/^ - e-^/^ (l - ^) , 
/L,(p) = e-^/^+'^^/^^e-^/^(l + ^) , 

/lib) = e-^/^-^^/^ ^ e-^/^ (l - , (27) 

where /ii and /i2 are chemical potentials related to the (approximately) conserved baryon 
and lepton numbers in the fourth generation. In expanding the exponents we have used 
the fact that baryon and lepton excesses are small. From ( pTj) we obtain 



gi=l s=l ' 
2 2 



gi=l s=l 



rf3 



EE/ 773^fMfuip)=L., (28) 



(27r 



which relate b and / to the density of excess baryons and leptons respectively. The sums 
are over the fiavor and spin indices in the fourth generation. 

One can now use (p^) and (^) in (^) and express products like fqj_{j)i)fqi{p2) as 
/^^(pi+P2)(l + 26) in decay and inverse decay terms. Subtracting the Boltzman's equation 
for the antiquarks from the equation for the quarks one gets the equation for the baryon 
number in the fourth generation: 

dB 1 

—^ + 3HB^ ^ -{nx -n''^){ri -T]){T{X ^ q^q^^q^h)) 
dt 2 
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^(3^4 + Li)n^{a{qiqi qJi)) 



6 riy 



2 (54 + L,){TiX ^ q,k)) + B,{T{X ^ q,q,)) 



- l{B,-Bs){T{q,^qsW)), (29) 

where {T{X — ^ ^4/4)) and (r(X q^qi)) are averaged over initial state degeneracies and 
summed over final state degeneracies, {ai^q^^q^ qih)) is summed over both initial and 
final state degeneracies and ~ is the photon number density. The correct sign for 
the term — n^(r7 — rj)(r(X q^q/^^qllA)) is obtained only after including the CP violating 
part of |M'(g4g4 ^ qih)? - Wiq^qi ^ ^4^4)1' H- 



The various terms in the r.h.s of ( PU[ ) are readily interpreted. The first term is the 
driving term for baryogenesis. It becomes small as nx — ^ ffx ^^'^ plays no role in freeze 
out. The second term comes from inverse decays and the third term comes from 2^2 
baryon number violating processes. These two terms can potentially cause a washout. 
The last term is the rate at which baryon number is drained out of the fourth generation 
into the third generation. We have ignored similar drainage terms to other generations 
because the mixing angles are smaller by one or two orders of magnitude. 

The 'washout terms' can be ineffective only if they are smaller than the Hubble dilution 
term SHB^. We must, therefore, have (using L4 ^ B4) 



1 cq 

3HB4 > -^B^{T{X ^ qiq^.Tdl)) ■ (30) 



"1 



"iHBii > B4ny{a{q4^q4^ qlU)) 
1 n] 
6 n. 

5.3 Range of parameters 

At the weak phase transition if ~ 3 x lO^^^Tc. Using our estimates of the decay widths 



and cross sections from Appendix B we can reduce the conditions (|30|) to 

M9\ 



|4 



lO"^*^ > 



10-16 > l9lf^^/2^-k (31) 



where k = 

J c 



When these inequalities are satisfied, there is no significant washout of the baryon 
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number and the net baryon to entropy ratio is 

Vb = TfTs 10 1^1 k^e. (32) 

We have taken (47raqcd)^ = 2 and = 100. Two parameters in ( P^D are bounded from 
above. The maximum value of e ~ 1 and the maximum value of rixcrit ~ Tj}. When these 
bounds are taken into consideration, the inequalities (^Tj) and (|I]) yield the following range 
of values for k and {gl"^: 



25 < A; < 80, 
10"^ < < 10-^•^ (33) 



Picking some value for k further constrains the range for and vice versa. Realistic 
values for rjs can be obtained with these values. For instance, taking k = 30, {gl"^ = 10^^ 
we obtain, 

r^B ~ 2.7 X 10-i°e. (34) 

For e close to 1 this falls within the range given by (|l|). Note that for k ~ 25, mx ~ 6.25 
TeV. This value is to be compared with the mass of the smallest string loops. Indeed if the 
mass per unit length of the strings is fi, a string loop of size i? ~ ^ has a mass of about 
where /3 is a numerical factor that takes into account the fact that loops are not 
exactly circular. Numerical simulations indicate that /5 ~ 9 |2^. If /i ~ (TeV) 2 then the 
mass of the smallest loops is about 36 TeV. Also note that r]B is insensitive to the number 
of Xs produced per string loop as long as NTj} > nxcrit- For = 10~^, /c = 30 we 
have nxcrit ~ lO^^T^^. 

The range of allowed values for the mixing angles is much wider. From (H) and (|33|) 
we can see that 

Now consider the decay of the fourth generation baryons and leptons. The decay widths 
\9w&i^^^\ Tc i^^Yieie gw is the weak gauge coupling). A lower bound on the mixing 
angles is obtained by requiring that the decays happen before nucleosynthesis. This means 
that the decay time should at most be Is. The corresponding limit on the mixing angles 
is: 9i4i3,9g4,3 > 10"^^ If we also require that ^ Qg^g,^ ^ lOO^^^j ~ lOOGg^gj, then 
we have 

10-'=^ < e,,,3 < 10-^ (36) 
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6 Conclusions 



We have shown that baryogenesis from the decay of hght scalar bosons is viable even 
at energies as low as the electroweak scale. This is interesting, since perturbative vio- 
lation of baryon number at low energies seems incompatible with the observed stability 
of the proton. However the minimal model of B, C and CP violation involving light 
scalar bosons can naturally have very small mixings between new heavy quarks with B 
violating interactions and the lighter quarks which shields the proton from B violating 
effects. Other phenomenological and cosmological constraints are shown to be satisfied. 
In particular the small mixings between the fourth family and the other families is shown 
to be sufficient for quarks in the fourth generation to decay into quarks of the lighter 
generations in a cosmologically acceptable time. 

Some members of the fourth family can be as light as 100 GeV. They can also be 
relatively long lived (decay time ~ 10~^s). It would be interesting to explore signatures 
of their existence in future experiments. In particular, if the lighter quark in the fourth 
family, the 6', happens to be lighter than the top quark then its dominant decay mode 
is the FCNC mode 6' — > fo + Z. One expects this decay mode to be explored in collider 
experiments of the near future. Signatures associated with new quarks with small mixing 
angles {Qq^q^ < lO""^) as predicted by this model, would seem to point strongly toward 
a mechanism of baryogenesis through perturbative B violation at the electroweak phase 
transition. 

The scalar bosons must be at least 25 times heavier than the electroweak scale. In our 
model they are produced copiously from the decay of loops of cosmic strings immediately 
after the electroweak phase transition. We show in Appendix A, it is possible to extend 
the standard model so that the smallness of the new mixing angles and the presence of 
the cosmic strings are justified simultaneously. 

Variations of this model can be conceived. The only necessary ingredients are topo- 
logical defects like cosmic strings and heavy baryons. The model has all the advantages of 
baryogenesis models where baryogenesis occurs after the electroweak phase transition in- 
cluding compatibility with the usual models of inflation. It is also viable as a baryogenesis 
model even if the electroweak phase transition is a second order phase transition. 
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Appendix A 

The model described in section 2 has an anomaly free global {B — L)^ symmetry in 
addition to the ususal B — L when the mixing angles between the fourth and lighter 
generation are put equal to zero and there is no coupling between the lighter generation 
quarks and leptons and the colored scalars Xi and X2. When the mixings and couplings 
mentioned above are non-zero but small, {B — L)4 is broken weakly. The smallness of 
these parameters is, therefore, not unnatural in the technical sense. Below, we describe a 
way of explicitly realizing this scenario as an effective low energy limit of a theory where 
the small {B — L)^ breaking terms come from operators of dimension 5 or higher and are 
suppressed by a large mass scale. 

We first define two C/(l) symmetries: 

4 3 
U = 5^(5 -L)i^B-L; V^{B-L),- J2{B - L)i . (37) 

i=l i=l 

The first is the usual B — L which we keep as a global unbroken symmetry of our theory. 
At a scale much higher than the electroweak scale one can conceive of a theory where V is 
a gauged symmetry. Consider for instance extending the model in section 3 by first taking 
away all terms violating V and then gauging V . (The terms involving the scalar X3 are 
not necessary for this discussion and can be discarded). In order that V be realized as 
a weakly broken symmetry in the low energy theory we can introduce a scalar field 
which is a singlet under all gauge symmetries except V . A vev for ^3 breaks V . 

The effective action for this theory can have a dimension 4 + n operator of the kind 
■^qihq^X'^ if the V charge of X3 is J^. M is a large mass suppressing this operator [h 
is the standard Higgs). When gets a vev, V is broken and we get the dimension 4 
mixing term { ^^^ YqAhqz- The smallness of the mixing results from suppression due to 
small coupling constants and factors of j- as well as the ratio 
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To actually get this operator we must introduce new fields and interactions that couple 
the fourth and the third generations. Since we are not interested in solving the hierarchy 
problem, scalars are cheap. We can introduce three more, X4, X5 and Xq with the 
following Yukawa couplings 

-^^Yukawa = fi'' (qiQaX^ + ^4/3X5 + ^3/4-^6) , (38) 

where qt and k denote quarks and leptons of the ith generation. We have chosen the 
same coupling g' for all the terms. The V charges of X4, X5 and Xg are and | and — | 
respectively. All are color triplets and SU{2)]y singlets. By suppressing the flavor and 
helicity indices we imply that all possible gauge invariant couplings are included in (|38|). 

Now suppose X4, X5 and Xg have masses of the order of M. Integrating them out 
one may obtain the dimension 6 operator q^q-dQ^h- This operator, a potentially dangerous 
candidate for proton decay, is not induced by renormalization at one loop and must be 
suppressed by at least a factor of ^^^sj^p where g is one of the couplings in (H). If we 
choose Igl'"^ < \g\'^ ~ 10~^, the proton decay problem is avoided for > 10^ GeV^. 

The mixing between the fourth and third generations can occur through the operators 
q^Dq^X^ (where D is the gauge-covariant derivative) or q^hq^X^. Figures 3a and 3b show 
typical leading order contributions to these operators. Clearly the V charge of X3 is |. The 
mixing is suppressed by the small number ■^^(■^^)^- For (X3) ^ 10^ GeV, M ^ 10^-^ 
GeV and |(y'||5'|' ~ 10~^, the mixing is 0^453 ~ 10"^*^. A similar mixing is obtained in the 
lepton sector. Mixings of this order are certainly small enough for the viability of our 
model. The mixings are also large enough for the baryon number in the fourth generation 
to be transported to the lighter generations in a cosmologically acceptable time. Indeed 
with this mixing the decay time of the fourth generation quarks and leptons is ~ 10~^s 
which corresponds to a temperature of about 1 GeV. By choosing (X3) ~ lO^GeV we also 
get the much needed cosmic strings at the electroweak scale as a bonus. 



Appendix B 

The dominant annihilation modes of X are XX GG, qq. Annihilations to lep- 
tons, Ws, Zs, Higgses and photons have much smaller rates because they are down by 
small coupling constants while the annihilation rates to quarks and gluons are enhanced 
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markedly by large color factors. We estimate the dominant annihilation processes in per- 
turbation theory. The lowest order Feynman diagram contributing to the annihilation to 
quarks is shown in figure 4a. There is a single gluon exchange. The squared amplitude, 
after summing over final state (spin, color and fiavor) degeneracies and averaging over 
initial state degeneracies, is 



|M(X ^qq)Y = — (47raqcd)'^sin^^, (39) 



where p is the 4 momentum of the X and 6 is the scattering angle in the cm frame. We 
have taken all the quarks to be massless. Feynman diagrams corresponding to annihilation 
to gluons are shown in figure 4b. The invariant squared amplitude is 



|M(XX^G'G')p = -y^^^i^ 



1 + cos» , (40) 



where is the angle between the two final state gluons in the rest frame of one of the 
incoming particles and p is the 4 momentum of the X in the cm frame. The large 
numerical factors in (|39|) and (^) come from color and fiavor sums. Note that we have 
four generations of quarks now. 

The thermally averaged cross section can be approximated by the zero temperature 
cross section with \p\ p^ Mx- We then obtain 

(a(XX^GG,gg))^%^. (41) 
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The XX annihilation rate is to be compared with the baryon number violating decay 
rates of Xs. The lowest order Feynman diagrams for these decays are shown in figure 
4c. The squared amplitude for the decay to two quarks (after summing over final state 
degeneracies and averaging over initial state degeneracies) is 



|M(X^g4g4)|'~4|(7pMi. (42) 

Once again we have taken the final state particles to be massless. The corresponding rate 
for a decay to an anti-quark and an anti-lepton is exactly the same (larger fiavor factor 
compensates for the smaller color factor): 

\M{X^7uU)\'^Mg\'M],. (43) 
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Approximating the thermally averaged decay width by the zero temperature decay 
width we get 

(r(X ^ qm,qiU)) ^ ^\9\'Mx • (44) 

Figure 4c shows the Feynman diagram corresponding to the leading order contribution 
to the process ^4^4 ^4/4. The invariant squared amplitude for a typical process is 

\M{q,q,^qr,)\'^9Q\g\'^, (45) 

where kg is the cm energy of a ^4 in the initial state. The thermally averaged cross 
section (summed over all initial and final state degeneracies) is approximated by a zero 
temperature cross section with /cq set equal to T^. The result is 

{a{qm^m)^^^j£^- (46) 
Finally the ^4 — >• q^W and I4 I3W decays (figures 4d) have the widths 

{Tiq,^q,W)) = ^QlJgw\'T,, 

{Tih^kW)) = j^Ql^\gw\'n, (47) 

where we have made the approximation that the final state particles are much lighter than 
the decaying particle and averaged over initial state degeneracies. We have also approx- 
imated the thermal averaging by taking all masses and momenta in the final expression 
to be of order Tc. Even with the limitations of the above approximations, the expressions 
in (H^) are useful as order of magnitude estimates of these decay rates. 
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Figure 2. Proton decay to pion and positron. 




Figure 3a. Dim. 6 operator q^Dq^XjX^ 
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Figure 3b. Dim. 6 operator hq^ ^3^3 
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Figure Ah.X annihilation to gluons. 
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Figure 4 a. X annihilation to quarks. 
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Figure 4 c. 2 to 2 baryon number violating process. 




Figure 4 d. Decays of q4 and I4 



